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Abstract 

The wave function of a spheroidal harmonic oscillator without spin-orbit interaction 
is expressed in terms of associated Laguerre and Hermite polynomials. The pairing gap 
. and Fermi energy are found by solving the BCS system of two equations. Analytical 

CNI ! relationships for the matrix elements of inertia are obtained function of the main quan- 

6J[)' turn numbers and potential derivative. They may be used to test complex computer 

codes one should develop in a realistic approach of the fission dynamics. The results 
given for the 240 Pu nucleus are compared with a hydrodynamical model. The impor- 
tance of taking into account the correction term due to the variation of the occupation 
number is stressed. 
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O 1 Introduction 

O 

By studying fission dynamics P one can estimate the value of the disintegration constant 
A of the exponential decay law expressing the variation in time of the number of decaying 
nuclei. The partial decay half-life T is given by T = rln2 = 0.693147/A. The potential 
energy surface in a multi-dimensional hyperspace of deformation parameters j3\, fa, fa 
gives the generalized forces acting on the nucleus. Information concerning how the system 
reacts to these forces is contained in a tensor of inertial coefficients, or the effective mass 
parameters {B^}. Unlike the potential energy E = E(/3) which depends on the nuclear 
shape, the kinetic energy is determined by the contribution of the shape change expressed 
by 

E^\±B,AP)ff (1) 

where Bij is the inertia tensor. In a phenomenological approach based on incompressible 
irrotational flow, the value of an effective mass B trr is usually close to the reduced mass 
/i = (AiA 2 /A)M in the exit channel of the binary system. Here M is the nucleon mass. One 
may use the Werner-Wheeler approximation 0. 

The microscopic (cranking) model introduced by Inglis jS] leads to much larger values of the 
inertia. By assuming the adiabatic approximation the shape variations are slower than the 
single-particle motion. According to the cranking model, after including the BCS pairing 
correlations |UEj, the inertia tensor is given by jEllZ|: 
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where H is the single-particle Hamiltonian allowing to determine the energy levels and 
the wave functions \u), u u , v v are the BCS occupation probabilities, E v is the quasiparticle 
energy, and Pij gives the contribution of the occupation number variation when the defor- 
mation is changed (terms including variation of the gap parameter, A, and Fermi energy, A, 



Similar to the shell correction energy, the total inertia is the sum of contributions given 
by protons and neutrons B = B p + B n . The denominator in equation (J2J) is minimum for 
the levels in the neighbourhood of the Fermi energy. A large value of inertia is the result 
of a large density of levels at the Fermi surface. As a result, in a similar way to the shell 
corrections, one can observe large fluctuations of Bu when the deformation or the number 
of particles are changed. 

In the present work we consider a single-particle model of a spheroidal harmonic oscillator 
without spin-orbit interaction for which the cranking approach allows to obtain analytical 
relationships of the nuclear inertia. Despite of the limited interest of this simple single- 
particle model, the result of the present work may be used to test complex computer codes 
developed in a realistic treatment of the fission dynamics based on the deformed two center 
shell model [Hj. The results illustrated for 240 Pu nucleus are compared with a hydrodynamical 
model. 

2 Nuclear shape parametrization 

The shape of a spheroid with semiaxes a, c (c is the semiaxis along the symmetry) expressed 
in units of the spherical radius Rq = tqA 1 ^ 3 may be determined by a single deformation 
coordinate which can be one of the following quantities: the semiaxes ratio c/a; the eccen- 
tricity e = (1 — a 2 /c 2 ); the deformation 5 = 1.5(c 2 — a 2 )/(2c 2 + a 2 ), or the quadrupolar 
deformation jHj e = 3(c — a) /(2c + a) which will be used in the following, and according to 
which the two oscillator frequencies are expressed as: 



and by taking into account the condition of the volume conservation io\io z = (wq) 3 where 
Twul = 41A" 1 / 3 MeV, one has 



dA/d/3i and dX/dpi): 





(3) 
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A particularly interesting value is e = 0.6 of a super deformed spheroid with the ratio c/a = 2. 



3 Spheroidal harmonic oscillator 

The eigenvalues £Q in units of Huq are given by 

E = hu L {n L + 1) + hu z (n z + 1/2) = hu [N + 3/2 + e{n ± - 2N/3)} (6) 

where the quantum numbers n± and n z are nonnegative integers. Their summation gives the 
main quantum number N = n± + n z . In units of hu Q one has a linear variation of the energy 
levels in function of deformation e. By including the variation of ujq with e, one obtains the 
analytical relationship 

d = Ei/hul = [N+ 3/2 + e(n± - 2iV/3)][l - e 2 (l/3 + 2e/27)]~ 1/3 (7) 

in units of TiuJq. Due to the Pauli principle, each energy level e*, with quantum numbers n± 
and N, can accomodate (7 = 2n± + 2 nucleons. One has a number of (N + 1)(N + 2) nucleons 
in a completely occupied shell characterized by the main quantum number N, and the total 
number of states for the lowest JV+1 shells is J2n=o( n + 1 )( n + 2 ) = (N+l)(N+2)(N+3)/3. 
For each value of there are N + 1 levels with n± = 0, 1, 2, ...N. When the deformation 
e > increases, a level with n± = decreases in energy and the one with n± = N increases. 
For some particular values of the deformation parameter there is a crossing of several levels 
in the same point leading to a degeneracy followed by an empty gap. If no spin-orbit coupling 
is considered for the vanishing deformation parameter, e — 0, one has the following sequence 
of magic numbers: 2, 8, 20, 40, 70, 112, 168, 240, ... If now e = 0.6 they become 2, 4, 10, 
16, 28, 40, 60, 80, 110, 140, 182, ... The known experimental values can be obtained only 
with a spin-orbit coupling included. 

The spin E contributes with positive or negative values (up or down) for every state with 
quantum numbers n z , n r = 0, 1, 2, ...n± and m = n± — 2n r , hence in a system of cylindrical 
coordinates (p, cp, z) the wave function [TDJ [TT] can be written as 

m = \n T mn z E) = C>)<M^n^)x(£) (8) 

Few examples of the quantum numbers n r mn z belonging to the lowest levels with N = 0,1,2 
are given in the table H 

The eigenfunctions are given by 

VC» = ^-N^e^L^in) = (9) 
= -^=e m * (10) 



2tt 



y 'a z y/a 

where are the associated (or generalized) Laguerre polynomials and H nz are the Hermite 
polynomials. The new dimension-less variables rj and £ are defined by 

"=5 (12) 
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Table 1: Quantum numbers of the lowest states of a spheroidal harmonic oscillator 



N 


n± 


n r 


n z = N — n± 


m = n_i_ — 2n r 


e for 5 = 
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£ = - (13) 

a z 



where the quantities 




A^jm a^J-^-^A^J^ (14) 



Mto i V u} i V Ma; 7 V w 



which depend on the nucleon mass, M, posses a dimension of a length. Their numeri- 
cal values, in fm, can be estimated by knowing that h 2 /M « 41.5 MeV-fm 2 and TiuJq = 
41A -1 / 3 MeV. The normalization constants 

(Nz? = , n /\ m (is) 

(n r + |m|)! 

(^) 2 = , (is) 

are obtained from the orthonormalization conditions 

«p)C>Wp = <W (17) 



i/} n '(z)i/; nz (z)dz = 6 n r nz (18) 



oo 
2tt 



$,* n ,(v?)$ m (^)^ = <Wm (19) 

One should take into account that the factorial 0! = 1! = 1. We shall substitute the wave 
functions il>™(r)) and ^ n -(0 i n the equation (J2J) of the nuclear inertia. 
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4 Nuclear inertia 



By ignoring the spin-orbit coupling the Hamiltonian of the harmonic spheroidal oscillator 
contains the kinetic energy and the potential energy term, V: 

T r/ * \ 1,,/ 22 2 2n 1* U .2 [(3 + e)r? + (3 - 2e)^ 2 l 

Now we are making some changes in the equation (J2J), first of all replacing the deformation 

(3 by £. 

One may assume jHl 13 El that only the leading term of the Hamiltonian, namely the 
potential written above, contributes essentially to the derivative, 

d JL * *L (2i) 

de de 

The contribution of Py, denoted by P £ for a system with one deformation coordinate, some- 
times assumed to be neglijible small, will be discussed in the last section. 

The derivative is written as 

= §[/*),+ AW?] (22) 

in which 

f = g ( g + 6 ) + 9 (2V 
Jl [27-e 2 (9 + 2e)] 4 / 3 1 ; 

f = 2 g ( 2g + 3 )~ 9 r24l 
72 [27-£ 2 (9 + 2e)] 4 /3 ^ 

For a single deformation parameter the inertia tensor becomes a scalar B £ whith a summation 
in eq. El performed for all states u, \x taken into consideration in the pairing interaction [T2*] . 



4.1 Pairing interaction 

We consider a set of doubly degenerate energy levels {e^} expressed in units of HuJq. Calcula- 
tions for neutrons are similar with those for protons, hence for the moment we shall consider 
only protons. In the absence of a pairing field, the first Z/2 levels are occupied, from a total 
number of n t levels available. Only few levels below (n) and above {n') the Fermi energy 
are contributing to the pairing correlations. Usually n! = n. If g s is the density of states at 
Fermi energy obtained from the shell correction calculation g s = dZ/de, expressed in number 
of levels per HuJq spacing, the level density is half of this quantity: g n = g s /2. 

We can choose as computing parameter, the cut-off energy (in units of TiuIq), Q ~ 1 3> A. 
Let us take the integer part of the following expression 

Qg s /2 = n = n' (25) 

When from calculation we get n > Z/2 we shall take n = Z/2 and similarly if n' > n t — Z/2 
we consider n' = n t — Z/2. 
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The gap parameter A = |G| J2k u kVk and the Fermi energy with pairing corellations A (both 
in units of hu^) are obtained as solutions of a nonlinear system of two BCS equations 

ri - n = V . €k ~ X (26) 
k % ^{e k - A) 2 + A 2 

2 kf 1 

G = £ 75 - A) 2 + A 2 (2T) 
where fcj = Z/2 — n + 1; kf = Z/2 + n'. 

The pairing interaction G is calculated from a continuous distribution of levels 

2 g(e)de , . 

(28) 



G JX-u ^(e- A) 2 + A 2 

where A is the Fermi energy deduced from the shell correction calculations and A is the gap 
parameter, obtained from a fit to experimental data, usually taken as A = 12/ y/Ahu)®. Both 
A p and A n decrease with increasing asymmetry (N — Z)/A. From the above integral we get 

2fif(A) In ( 2 4) (29) 



G av ' VA 
Real positive solutions of BCS equations are allowed if 

f?^V : (30) 

i.e. for a pairing force (G-parameter) large enough at a given distribution of levels. The 
system can be solved numerically by Newton-Raphson method refining an initial guess 

A = {n s e d + n d e s )/(n s + n d ) + G(n s - n d )/2 

A 2 = n s n d G 2 - (e d - e s )/4 (31) 

where e s , n s are the energy and degeneracy of the last occupied level and e d , n d are the same 
quantities for the next level. Solutions around magic numbers, when A — > 0, have been 
derived by Kumar et al. |14j . 

As a consequence of the pairing correlation, the levels situated below the Fermi energy are 
only partially filled, while those above the Fermi energy are partially empty; there is a given 
probability for each level to be occupied by a quasiparticle 



1 _ e fc - A 



(e k - A) 2 + A 2 



(32) 



or a hole 

< = 1 - v 2 k (33) 

Only the levels in the near vicinity of the Fermi energy (in a range of the order of A around it) 
are influenced by the pairing correlations. For this reason, it is sufficient for the value of the 
cut-off parameter to exceed a given limit Q 3> A, the value in itself having no significance. 
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4.2 Variation with deformation 

The following relationship allows to calculate the effective mass in units of h 2 / (hu^) 







h 2 2 ^ 



(34) 



Matrix elements are calculated by performing some integrals 

-2 



where 



n' z n' r m'\h{e)r) + f 2 {e)e\n z n r m) = S m r m N^ N™ N n > z N n _ ■ [ 
fi / rfW m|+1 e- r '47' mt 1 fo) / #e-* ^ (O^n. (0 + 

JO J — oo 

^0 J —oo 

J —oo 

d vv ^L^( v )L^( v )=5 Knr 



(35) 
(36) 



so that 



(n' g n' r m'\h(e)ri + f 2 (e)£ 2 \n z n r m) = 6 mlm N^N^N n ,N nz - 

JO J —oo 

Next we can use the relationships |16j : 

(n r + |m|)! 



d^^+VLjJ'^LC'^) = <W2n r + H + 1)- 



(37) 



/oo 
-oo 



: ( n z + n) + <5n'n z +2(n z + l)(ri 2 + 2) + (5 



(38) 

which were obtained by using the recurrence relations and orthonormalization conditions for Her- 
mite polynomials and associated Laguerre polynomials. 

L*(x) = L* +1 {x) - L k +\(x) ; ( n + l)L* +1 (x) = [(2n + fc + l)-x]L*(x)-(n + fc)L*_ 1 (x) (39) 

H n+1 (x) = 2xH n (x) - 2nH n ^{x) (40) 
with particular values Lq(x) = 1, L\(x) = 1 + k — x, Hq{x) = 1, H\{x) = 2x. 

Eventually from the sum of equation (|34jl one is left with an important diagonal contribution and 
two nondiagonal terms 



B £ l = ~;5n' r n r 5 m 'm iM 2n r + H + X ) + h(n z + 1/2)]' 



(UuVu) 



u=ki 



-D £ 2 = ~:0 n ' r n r 0m'm ~7r\ n z + l)(n 2 + 2)—^ — , F vs n > z 



+2 



(41) 
(42) 
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Hn 9 e X ST fit 1\ + X //IQN 

where fcj and fc/ have been defined above. In order to perform the summations of the nondiagonal 
terms for a state with a certain v (specified quantum numbers n z n r m) one has to consider only the 
states with /i ^ v and n' r = n r ; m' = m for which = n 2 + 2 or ra^ = n z — 2 respectively. Finally 
one arrives at the nuclear inertia in units of h 2 /MeV by adding the three terms and dividing by 



4.3 Hydrodynamical formulae 

There are several hydrodynamical formulae ^2] of the mass parameters. For a spherical liquid drop 
with a radius #0 = 1.2249 A 1 / 3 fm one has 

2 2 f> 2 ft 2 

B irr (0) = —MAR 2 , = ^1.2249 2 ^ 5 / 3 fm 2 = 0.0048205^ 5/3 (44) 

1 ; 15 15 41.5MeV-fm 2 MeV V ; 

When the spheroidal deformation is switched on it becomes 

^w-^ V.^w-^ (45) 

Good results for the fission halflives of actinides have been obtained by using an inertia larger by 
about an order of magnitude. 

One can also employ a formula with a fit parameter k 

BP h (r)=v + k(B; rr (r)-v) (46) 

where k shows how much deviates the flow of nuclear matter from an irrotational one. By substi- 
tuting the above expression of the irrotational term, one obtains 



1 + t lZexp' 



15 r V d 

with parameter values of d = do = Rq/2A52 and k = 11.5, or d = 2d$ and k = 6.5. 



(47) 



5 Results 

The main result of the present study is represented by the equations Q41H43JI . which could be used to 
test complex computer codes developed for realistic single-particle levels, for which it is not possible 
to obtain analytical relationships. Nuclear inertia of 240 Pu calculated with the equation (|44[) for a 
spherical liquid drop and with eq. 1)45 Jl for spheroidal shapes is illustrated in figure ^ One can see 
how B irr (0) increases when the mass number of the nucleus is increased. The irrotational value 
B l "'(e) monotonously increases with the spheroidal deformation parameter e. Due to the fact that 
in this single center model the nucleus only became longer without developing a neck and never 
arriving to a scission configuration when the deformation is increased, the reduced mass is not 
reached as it should be in a two center model [2]. 

The cranking inertia of the spheroidal harmonic oscillator calculated by using the analytical rela- 
tionships (|41H43|) and the correction given in the next section shows very pronounced fluctuations 
which are correlated to the shell corrections (calculated with the macroscopic-microscopic method 
[T3*] ) plotted at the bottom of the figure ^ 
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Figure 1: Top: comparison of the effective mass (in units of h 2 /MeV) calculated by using the 
cranking model for the proton plus neutron level schemes, only for neutrons, as well as for the 
irrotational spheroidal and spherical shapes of 240 Pu. Bottom: shell corrections for neutrons 
and protons, only for neutrons, pairing corrections, and shell plus pairing corrections. 

5.1 Variation of the gap parameter and Fermi energy with defor- 
mation 

Now we can calculate the correction term as 



2h 2 1 



A— ] 

de , 



dA 



dX dA 



( £ „-A)^_j + 2Afo-A)-_ 

-2A 2< ^{v\dV/de\u) 
de 

-2A(e v - X)^—(u\dV/de\u) 
de 

In figure |21 we plotted the variation with deformation of the solutions of BCS equations for Fermi 
energy A (bottom) and the gap parameter A (top) of the proton and neutron level schemes for 
240 Pu nucleus. The dotted line at the value 0.117 corresponds to A. Their derivatives with respect 
to the deformation parameter are given in figure |31 For superdeformed nuclei with e > 0.5 the 
oscilllation amplitudes of d\ n /de are approaching their maximum values of about 2 units. In the 
same range of the deformations the inertia is also larger as a result of the increased density of levels 
at the Fermi surface. 

The result displayed in figure |I] shows the important contribution of the neutron level scheme, 
P en (dotted line), reflecting the larger density of states at the Fermi energy, compared to the 
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Figure 2: The variation with deformation of the solutions of BCS equations for Fermi energy 
A (bottom) and the gap parameter A (top) of the proton and neutron level schemes for 240 Pu 
nucleus. The energies are expressed in units of TiuJq = 6.597 MeV. The dotted line in the 
upper part corresponds to A = 0.117. 

proton term P ep (dashed line). Their sum is a positive quantity, contributing to an increase of 
the nuclear inertia. In a dynamical investigation using the quasiclassical WKB approximation, 
the quantum tunnelling penetrability depends exponentially on the action integral, in which the 
integral contains a square root of the product of mass parameter and deformation energy. This 
exponential dependence amplifies very much any variation of the inertia. Consequently, the term 
Pij should be considered in calculations. A similar conclusion was drawn from a study of a realistic 
two-center shell model |17j . 

6 Conclusions 

By using the wave functions of the spheroidal harmonic oscillator (the simplest variant of the Nilsson 
model) without spin-orbit coupling one can obtain analytical relationships for the cranking inertia. 
Consequently the result may be conveniently used to test complex computer codes developed for 
realistic two center shell models. 

Unfortunately this single center oscillator is not able to describe fission processes reaching the 
scission configuration or ground states with necked in or diamond shapes. When the deformation 
parameter is increased the nucleus became longer and longer without developing a neck and reaching 
the touching point configuration. In this way it is not possible to obtain at the limit a nuclear 
inertia equal to the reduced mass of the final fragments in a process of fission, alpha decay or 
cluster radioactivity. 

As expected, in agreement with the results obtained by other authors, the cranking inertia is larger 
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Figure 3: The derivatives with respect to deformation of the solutions of BCS equations for 
Fermi energy A and the gap parameter A of the proton and neutron level schemes for 240 Pu 
nucleus. The energies are expressed in units of TiuJq = 6.597 MeV. 




Figure 4: Contribution, P e , to the mass parameter of the occupation number variation with 
deformation for 240 Pu nucleus expressed in units of h 2 /MeV. 
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than the hydro dynamical one for a spheroidal shape which is higher than that of a spherical nucleus. 

Acknowledgment This work was partly supported by Deutsche Forschungsgemeinschaft, Bonn, 
and by the Ministry of Education and Research, Bucharest. We acknowledge also the support by 
Prof. S. Hofmann, Gesellschaft fur Schwerionenforschung (GSI), Darmstadt. 

References 

[1] D. N. Poenaru and W. Greiner, in Nuclear Decay Modes, (Institute of Physics Publishing, 
Bristol, 1996) p. 275. 

[2] R. A. Gherghescu, D. N. Poenaru and W. Greiner, Phys. Rev. C 52 (1995) 2636; Z. Phys. A 
354 (1996) 367. 

[3] D. R. Inglis, Phys. Rev. 96 (1954) 1059; 97 (1955) 701; 103 (1956) 1786. 
[4] J. Bardeen, L. Cooper and J. Schrieffer, Phys. Rev 108 (1957) 1175. 
[5] S. T. Belyaev, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 33 (1959) No. 2. 
[6] D. R. Bes, K. Dan. Vidensk. Selsk. Mat. Fys. Medd. 31 (1961) No. 11. 

[7] M. Brack, J. Damgaard, A. Jensen, H. C. Pauli, V. M. Strutinsky and G. Y. Wong, Rev. Mod. 
Phys. 44 (1972) 320. 

[8] R. A. Gherghescu, Phys. Rev. 67 (2003) 014309. 

[9] S. G. Nilsson, Dan. Mat. Fys. Medd. 29 (1955) No. 16, 1. 

[10] J. Damgaard, H. C. Pauli, V. V. Pashkevich and V. M. Strutinski, Nucl. Phys. A 135 (1969) 
432. 

[11] D. Vautherin, Phys. Rev C 7 (1973) 296. 

[12] M. Bolsterli, E. Fiset, J. R. Nix and J. Norton, Phys. Rev C 5 (1972) 1050. 
[13] V. M. Strutinsky, Nucl. Phys. A 96 (1967) 420. 

[14] K. Kumar, B. Remaud, P. Aguer, J. Vaagen, A. Rester, R. Foucher and J. H. Hamilton, Phys. 
Rev C 16 (1977) 1235. 

[15] A. Sobiczewski, E.Ch.A.Ya. 10 (1979) 1170. 

[16] E. W. Weisstein, From MathWorld A Wolfram Web Resource 

http: / /mathworld. wolfram.com/LaguerrePolynomial. html , 
http: / / mathworld.wolfram.com/HermitePolynomial.html| 

[17] W. Schneider, J. A. Maruhn, W. Greiner, Z. Phys. A 323 (1986) 111. 



12 



